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1. Introduction 

Varying temperature and baryon density, it is possible to realize several different phases of 
nuclear matter. A qualitative phase portrait of admissible phases, on temperature-baryon density 
plane, is shown in Fig. 1, where po = 0.167 fm 3 is the normal baryon density (e.g., [1]). The 
variation of temperature and baryon density can be achieved in hadron-hadron and heavy-ion col¬ 
lisions. The corresponding values of these thermodynamic variables could also exist at the early 
stages of the Universe evolution or in the cores of neutron stars. There can occur phase transitions 
of first order, second order, as well as crossovers [2-4]. 


T(MeV) 



Figure 1: A qualitative phase portrait of admissible phases, on temperature-baryon density plane. 


Phase transitions are known to be difficult for description because of the absence of small 
parameters in the transition region. This especially concerns the description of phase transitions 
in microscopic theory, where calculations are possible only by resorting to a kind of perturbation 
theory. However, perturbation theory usually results in divergent series that, in the best case, have 
meaning in the limit of asymptotically small parameters, while the physical parameters could be 
rather large, or even infinite. The notorious question is: How it would be possible to extract in¬ 
formation from perturbative series, derived for asymptotically small parameters, for the values of 
finite and large parameters? 

There exist methods for effective summation of divergent series, such as Pade summation [5] 
and Borel summation [6]. However, the former exhibits a number of deficiencies, including the 
appearance of spurious poles, while the second is quite complicated and requiring the knowledge 
of large-order terms in power law expansions. Both these methods not always are applicable, as is 
discussed in [7]. 

In the present report, we describe an original approach to treating divergent perturbative series, 
extracting from them meaningful answers, providing good accuracy, with being much simpler and 
more general than Pade-Borel summation. We illustrate the approach by the example of a phase 
transition in multicomponent field theory at finite temperature. 
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2. Optimized perturbation theory 

The first step of the approach is optimized perturbation theory, based on the definition of 
control functions reorganizing divergent series to convergent ones. The optimized perturbation 
theory was advanced in 1973, in Thesis [8], submitted for publication in 1974, and published 
[9,10] in 1976. The power of this method was illustrated by treating some anharmonic models 
and strongly anharmonic quantum crystals [8-18]. Later this theory has been applied to numerous 
models, under different guises and under different names, such as modified perturbation theory, 
variational perturbation theory, renormalized perturbation theory, oscillator representation, delta 
expansion, optimized expansion, nonperturbative expansion, and so on (e.g., [19-24]). All these 
works have used the variants of the same idea [8-18] of introducing control functions renormalizing 
divergent perturbative series into convergent series. In this section, we briefly delineate the idea of 
optimized perturbation theory, as advanced in [8-18] and reviewed in [25,26]. 

Suppose we are interested in finding a function f(x ) satisfying a complicated equation that 
cannot be solved exactly, but can be treated only by a kind of perturbation theory. To explain the 
main idea, we consider here, for simplicity, a real function of a real variable. A generalization to 
complex functions and variables is straightforward. 

Let perturbation theory give a divergent perturbative sequence {/^(jc)}, with k = 0. 1.2,... 
being an approximation-order index. The basic idea of optimized perturbation theory is to intro¬ 
duce a set {uk{x)} of control functions that would allow us to reorganize the divergent perturbative 
sequence {fk(x)} into a convergent sequence {Ffx.uf}, where = u/fix). A sequence is conver¬ 
gent if and only if it satisfies the Cauchy criterion: For each positive e, there exists an order k e , 
such that 

| F k+p (x,u k+p )-F k (x,u k )\ < £ , (2.1) 

for k > k £ and any positive p = 0,1,2,.... 

Control functions can be introduced in several ways, three of which are the most common: (i) 
through initial conditions, (ii) through variable changes, and (iii) through functional transforma¬ 
tions. 

Introduction of control functions through initial conditions 

The simplest illustration of this way is when the sought function f{x) is a solution of a func¬ 
tional equation 

£\f{x)\ = 0 • (2.2) 

Then, starting with an initial approximation Fq(x,u), including a control function u, it is admissible 
to represent the functional equation ( |2.2| ) as an iterative procedure 

F k (x,u k ) = F k _x{x,u k -\) + S\F k ^\{x,u k ^j] . (2.3) 

As a rule, physical systems are characterized by their Hamiltonians, or Lagrangians. Say, H 
is a Hamiltonian of a complicated system that can be treated only by means of perturbation theory. 
Taking for the zero approximation a simple Hamiltonian Hq{u), containing control functions, one 
can rewrite the system Hamiltonian as 

H = H q (u) + [H-H 0 (u)\. (2.4) 
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Then, by perturbation theory with respect to the difference (H — Ho), one obtains higher approx¬ 
imations, depending on the considered problem, either for wave functions or for Green functions. 
Knowing the latter, one can calculate the corresponding approximations for observables A{x) as 
the averages 

F k (x,u k ) = (A(x)) k , (2.5) 

defined for the related k-order approximate wave, or Green, functions. 

Introduction of control functions through variable changes 
It is possible to change the variable x through the relations 

x = x k (z, Uk) , z = Zk(x, Uk) , (2.6) 

involving control functions, thus, getting 

fk{x) = fk{x k (z,u k )) . (2.7) 

The latter expression can be expanded in powers of the new variable z, so that 

fk(xic(z,u k ))~f k (z,u k ) (z —>■ 0) . (2.8) 

With the inverse variable change ( |2.6[ ), one has 

fk(z,u k ) =f k (z k (x,u k ),u k ) . (2.9) 

This allows us to define 

F k (x, u k ) = f k (z k (x, u k ), u k ) . (2.10) 

In troduction of control functions through functional transformations 

The sought function can be subject to a transformation containing control functions, 

f(u)f(x)=F(x,u) , (2.11) 

with the inverse transformation 

f(x) = T~ x (u)F (x,m) . (2-12) 

Then we define 

F k (x,u k ) = T(u k )f k (x) . (2.13) 

Formulation of equations for optimal control functions 

After control functions are incorporated into F k (x. u k ), it is necessary to formulate explicit 
equations for their calculations. By their meaning, the control functions are to be defined in such 
a way that to induce convergence for the sequence {F k (x,u k )}. Since convergence is characterized 
by the Cauchy criterion (fPl), the optimal control functions, in the spirit of optimal control theory, 
can be defined as the minimizers of the Cauchy cost functional 

%[u\ = \ Y,\ F k+ P (x,u k+P ) -F k (x,u k )\ 2 , (2.14) 

2 * 
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whose mi ni miz ation provides the fastest convergence of the sequence {F k (x, u k )}. That is, we need 
to look for the minimal value 


min| F k+p (x,u k+p )-F k (x,u k )\ 


(2.15) 


of the difference F k Vp — F k for any given p. 


The minimization condition ( |2.15| ) involves two control functions u k+p and u k , which makes 
it impossible to define both of them simultaneously. Hence, we cannot find the exact absolute 
minimum of the Cauchy cost functional ( |2. 14 ), but we can try to find its approximate minimum. 
Assuming that the control functions u k+p and u k , and the related terms F k p and F k are close to 
each other, we can express Fk +P as 


Fk+ P (Fi tt-k+p) ~ Fk+ P (Fi Mit) T dll (fik+P ^k) 


Then minimization (2.15) becomes 


min 

U 


F k+p (x,u k )-F k (x,u k ) + 


dF k (x,u k ) 

du k 


(fik-\- P ti k J 


(2.16) 


(2.17) 


Depending on the relation between the difference F k+p — F k and the term containing the deriva¬ 
tive, there can be two cases. When the derivative term is smaller than the difference term, then 


minimization (2.17) is approximately satisfied under the minimal difference condition 


F k+p (x, u k ) - F k {x , u k ) = 0 . 


(2.18) 


But when the difference term is smaller than the derivative term, then minimization ( |2.17| ) is ap¬ 
proximately valid under the condition 


<*»-*) ^ = 0 


(2.19) 


This has to be understood as the minimal derivative condition 


dF k (x.u k ) 

du k 


( 2 . 20 ) 


provided the latter possesses a solution. In case there are no solutions, one has to set u k+p = u k . 

As is evident, the minimal difference and minimal derivative conditions are absolutely equiva¬ 
lent. Of course, in particular cases, one of them can yield a better accuracy than the other. However, 
in general, it is impossible to conclude that one is preferable to the other. 

In this way, the optimized perturbation theory [8-18] consists of the following steps. The di¬ 
vergent perturbative sequence {f k (x)} is reorganized into the sequence {F k (x.u k )} incorporating 
control functions u k = u k (x) making the latter sequence convergent. Control functions can be in¬ 
troduced in three ways, through initial conditions, through variable changes, or through functional 
transformations. Explicit equations for the control functions are derived from the minimization of 
the Cauchy cost functional. Approximate minimization can be done by means of either minimal 
difference or minimal derivative conditions. 
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3. Self-similar approximation theory 

Optimized perturbation theory has been used for various physical systems. It is not our aim 
here to give a review if these numerous applications. Just let us mention a couple of review-type 
articles [25,26], where further citations can be found. Despite a variety of very successful and wide 
applications of optimized perturbation theory, several questions remained unanswered: 

(i) How it would be possible to improve accuracy within the given number of perturbative 
terms? 

(ii) What is a necessary condition that the Cauchy cost functional could reach its absolute 
minimum, that is zero? 

(iii) How to decide which of the ways of introducing control functions would be the best one, 
when there are several such ways, say, by choosing different initial approximations? 

(iv) Could it be feasible to check the stability of the calculational procedure, when no exact 
solutions are available, that would allow for the explicit comparison of these exact solutions with 
the obtained approximations? 

(v) Is it possible to define general approximants, enjoying a fixed prescribed structure, extrap¬ 
olating the series, derived for an asymptotically small variable x —>• 0, to the arbitrary values of this 
variable from the whole interval [0,°°)? 

All these questions have been answered in self-similar approximation theory advanced in [27- 
33]. The main idea of this approach is to reformulate perturbation theory to the language of dy¬ 
namical theory, considering the approximation order k as discrete time, so that the approximation 
sequence {F k (x, u k )} be isomorphic to the trajectory of a cascade. Then the effective sequence limit 
will correspond to the cascade fixed point. And the control of the calculational procedure stability 
will be equivalent to the analysis of the dynamical system stability. 

Let us define the expansion function x = x k (<p) by the reonomic constraint 


Fo(x,u k (x)) = (p , x = xk(q>). 


(3.1) 


Introduce the endomorphism 

y k {q>) = F k (x k (q>),u k (x k (q>))) 
that, owing to constraint ( |3.1| ), enjoys the initial condition 

yoiv) = (p ■ 


(3.2) 


(3.3) 


The inverse to this endomorphism is 

F k (x,u k (x)) =y k (F 0 (x,u k (x ))) . 


(3.4) 


In terms of this endomorphism, the Cauchy cost functional (]2. 14) takes the form 


^ k 

This functional is exactly zero, provided that 

yk+ P (<p) =y/c((p) 


^p\u\ = \ LlM+p(p)-Tf(<P )| 2 • 

^ b 


(3.5) 


(3.6) 
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for all k > 0. In particular, for k = 0, we have 


y P M = yoM = (p ■ 


Combining (3Y) and (|T7|) yields the functional self-similarity relation 


(3.7) 


y k+P {(p) =yk(y P {(p )) • 


(3.8) 


possesses the semi-group property y k ° y p =y k+p , it can also be called 


Since transformation 
group self-similarity. 

Thus, the self-similarity relation (3.8) is a necessary condition for the Cauchy cost functional 
to be zero. Although it is not a sufficient condition. The family of the endomorphisms {y k }, with 
the group relation ( |3.8| ), forms a dynamical system in discrete time, termed cascade. By con¬ 
struction, the cascade trajectory {yk{tp) : k = 0,1,2,...} is bijective to the approximation sequence 
{F k (x,u k (x)) :k = 0,1,2,...}. 

The bijectivity of the cascade trajectory and approximation sequence means the following. If 
there exists the limit 


where 


then there also exists the limit 


such that 


y*((p(x)) = jimy k ((p) , 

k-> °° 

(3.9) 

(p(x) = lim F 0 (x,u k (x)) , 
k— 

(3.10) 

F*(x) = lim F k (x,u k {x)) , 
k— 

(3.11) 

F*(x)=y*(cp(x)). 

(3.12) 


Note that the existence of a trajectory limiting point y* guarantees the existence of a sequence 
limit F*, however this does not guarantee that F* necessarily corresponds to the sought function 
fix). Such a correspondence is an assumption typical of calculational procedures dealing with non¬ 
linear problems [34,35], for which the accuracy of approximations at each step cannot be explicitly 
established. 

For a dynamical system, the existence of a limiting trajectory point is equivalent to the exis¬ 
tence of a stable fixed point. Therefore, if the cascade trajectory {y&}, with increasing k, tends to a 
limiting point y*, the latter is a fixed point, such that 


yk(y*)=y*. (3.13) 

It is more convenient to deal with a dynamical system in continuous time, instead of a system in 
discrete time. This can be realized by embedding the approximation cascade into an approximation 
flow, 

{yk((p) '■ k € Z + } C {yr(<p) : 1 F ®+}> (3-14) 

with the flow trajectory passing through all points of the cascade trajectory: 

3 >t{<p)=yk(<p) ( t = k ). (3.15) 
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For the dynamical system in continuous time, it is straightforward to write down the flow 
evolution equation that is the Lie equation 


dt 


yt(<p) = v(y t (<p)) 


(3.16) 


where the right-hand side is the flow velocity 

v((p) = 


d t s 

Tt y,w 


J t =0 


Integrating ( 3.16 ) between a given point of the cascade trajectory y k = y k ((p) and a point 
y*k = yl(<p). we get the evolution integral 


f 

■hk 


y*k dy 


n{y) 


= 4, 


(3.17) 


in which t k is the time of motion from y k to y* k , while v k is the flow velocity on this time interval. 

Remembering that the cascade is embedded into the flow, the flow velocity, near the time t = lc, 
employing the Euler discretization, can be expressed through the cascade velocity 


v k {(p) = F k+ i(x k ,u k )-F k {x k ,u k ) + (u k+ 1 - u k ) tt— F k (x k ,u k ) , 

dujc 


(3.18) 


where x k = x k ((p) and u k = u k {x k ). Invoking the bijective relation between y k and F k , the evolution 
integral can be represented as 

r F k dxp 

Jf ,* Vfc(qj) 


= 4 • 


(3.19) 


The motion time t k can be treated as an additional control function. In the simplest cases, it can be 
set to one or l/k or defined through additional conditions [25,26]. 

If v k were zero, then y* k would be an exact fixed point of the cascade. Unfortunately, it is 
difficult to set v k zero, since it contains two, yet unknown, control functions. But we can require 
the minimal possible velocity, in that way defining the control functions u k = u k {x) by the condition 


min 

Uk 


F k+ \(x, u k ) - F k (x, u k ) + [u k+ \ - u k ) F k (x, u k ) 


(3.20) 


This minimization is equivalent to minimization (2.17), with p = 1. Analogously to the previous 
consideration, an approximate minimization can be done by one of the conditions (2.18) or (2. IS). 
Condition (2.IS) seems to be more convenient, which results in the velocity 

v k ((p) = F k+ i(x k ,u k ) -F k (x k ,u k ) . 

Employing this in the evolution integral ( 3.19| ) yields the renormalized approximant F k . 

Defining the control functions from the minimization of the cascade velocity implies that y k ((p) 
is an approximate fixed point, or quasi-fixed point. Respectively, F k (x,u k ) corresponds to an ef¬ 
fective limit f k (x) that is named the self-similar approximation of fix). In this way, we have the 
correspondence 

yt(<p( x )) = F k (x,u k (x)) -A fl (x) . (3.21) 
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The improvement of the accuracy of the self-similar approximation f£(x), as compared to 
the optimized approximation Ffix, lit (x)), is due to the following reason. The minimization of 
the Cauchy cost functional is equivalent to the minimization of the cascade velocity, which gives 
the optimized approximant Ffix,ufix)). However, the cascade velocity is not exactly zero. The 
Lie equation ( 3.1 6| ) describes the motion from the given optimized approximant Ffix,ufix)) to the 
quasi-fixed point F£(x,ufi) that improves the accuracy of the former approximant. 

As is mentioned above, to represent the effective limit of the approximation sequence, the fixed 
point has to be stable. The stability of the procedure here coincides with the stability of motion of 
the dynamical system, which is characterized by the map multiplier 


vfi<p) = -^yfi<p) ■ 

The motion at the point yfitp ) is locally stable, provided that 

\vfiyfi<p))\ < i • 

The multiplier at the quasi-fixed point is 

/ 4 *(<P) = Ffiyfiv)) ■ 

The quasi-fixed point is stable, when 

1 / 4 ( F C(x, u fix))\ < 1 , 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


where relation ( |3.21| ) is used. 

Because in the treated case, the fixed point is a function of x, it is possible to consider the 
maximal multiplier 

A4* = SU P (<P)I = sup jjufiF k *(x,ufix)))j . (3.26) 


Then we say that a quasi-fixed point is uniformly stable, when 

1/4*1 <1- 


(3.27) 


The analysis of the procedure stability makes it possible to answer the question on which of 
the procedures is preferable, when there are several admissible procedures differing by the way of 
introducing control functions. For example, it is possible to introduce control functions by different 
initial approximations, as has been analyzed for anharmonic models [25]. For strongly anharmonic 
quantum crystals, it is possible to choose different initial approximations, say, Hartree or Hartree- 
Fock [17,36,37]. Or one can introduce control functions by different changes of variables [6], 

The answer is: That procedure is preferable that is more stable, since a more stable procedure 
is assumed to be faster convergent [25,26]. 

Finally, we give the answer to the problem whether it is feasible to construct general expres¬ 
sions extrapolating the series in powers of an asymptotically small variable x —>• 0 to its arbitrary 
values in the whole range x € [0,°°). 

Suppose, the sought function can be found only for an asymptotically small variable, 


f(x)-ffix) (*->0), 


(3.28) 
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where it is given by the asymptotic expansion 



Such series are usually divergent for any finite x. 

It is convenient to consider the normalized function defined by the ratio 


fk(x ) 
fo(x) 



(3.29) 


(3.30) 


that, by construction, satisfies the limit 


fkix) 
^0 fo(x) 


= 1 


Control functions can be introduced by applying to series (3.30) the method of fractal trans¬ 
forms [26,38-45] defined as 

FAx,u) = x u 
^ fo(x) 


Then, following the self-similar approximation theory by accomplishing several times the renor¬ 
malization procedure, we come, depending on the available boundary conditions, to one of the 
following approximants. 


Right root approximants 


m =((... (1 +Aijc)"‘ +A 2 x 2 Y 2 + ... +A k x k J k (3.31) 

that can be used, when a number of terms in the large-variable expansion x —)■ °° are known. All 
parameters Aj and n-, are uniquely defined through this expansion [26,41,43,44]. 

Left root approximants 


fo(x) 


^(- ■ • (1 f-Aix)" +A 2 X 2 ) + ■ • • + Akxl^j , (3.32) 


in which the sole power n* is defined from the large-variable behavior x while all parameters 
Aj are found from the accuracy-through-order procedure after re-expanding ( [3.32 ) in powers of 
x —> 0 and comparing this with the initial expansion ( |3,29| ). We may note that ( 3.32 ) is a particular 
case of (|3.31|), with 


7+1 


(7 = 1 , 2 , 1 ) , 


which involves all Aj in the definition of the large-variable amplitude [7,46]. 
Continued root approximants 


fkW 

fo{x) 


(1 + A\x (1 A 2 x ... (1 -\-Akx) s y ...) , 


(3.33) 
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in which the power s is prescribed by the large-variable behavior, while all parameters A,- are given 
by the accuracy-through-order procedure at x —> 0. In the particular case of s = — 1, approximants 
( |3.33 ) are reduced to continued fractions and, hence, to Pade approximants [47]. 


Exponential approximants 

A'W 

/oM 


= exp(b\x expfax.. .exp(b k x ))...) 


(3.34) 


where the control functions b, are defined by additional conditions [26,40-42], like the minimal 


difference condition ( |2.18[ ). More elaborated variants of defining these control functions are also 
possible [26], e.g., 

a„(l + a\) 


bn — /, . 2 '\ 

na n -\(l + affi 

where a n are the parameters of expansion ( |3.29| ). 
Factor approximants 


(n = 1,2,...,*) , 


in which 


fk(x) 

fo(x) 


N k 

n( 


1=1 


N k 


k/2, k = 2,4,... 
(* +1)/2, * = 3,5,... 


(3.35) 


and all parameters A, and /?, are defined from the re-expansion procedure at x —> 0, equating the 
like-order terms [48-52]. 


4. Multicomponent field theory 


Let us illustrate the application of self-similar approximation theory to describing a phase 
transition in IV-component <p 4 field theory in d-dimensional space. The Hamiltonian of this field 
theory is 


H[cp} = 



d(jo(x) 

l 2 

dx 


2 m 2 j. . A a. , I , 

+ -y <P~(x) + y W } dx , 


where the standard notations are employed: 

<p(x) = {< p„(x ) : n = 1,2,... ,7V} , x = {x a : a = 1,2 


‘ P 2 (x ) = £ %(x) 

n= 1 


d(p(x) 

dx 


f f \ d%(x) 

h hi l dx « 


Hamiltonian ( p~I| ) is invariant under the reflection 

%(x) -%(x) (n = 1,2,...TV) , 


(4.1) 


(4.2) 


so that H[—(p] = H[<p\. This means that the statistical average of the field is zero: (<p) = 0. The 
thermodynamic potential 

F[((p)} = -7TnTr<T™ , (4.3) 
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where T is temperature and /3 = 1/7", is also invariant under reflection (|4~2|). 

It turns out that there exists a critical temperature T c , such that above this temperature, the 
order parameter (tp) is zero, while below T c it is nonzero: 

(<?(*))= 0 (T > T c ), 

0 (T < T c ), 

which implies that 

F[(cp(x))^0]<F[((p(x)) = 0], 

for T < T c . The zero {(p) means that all (%) are zero. And nonzero ((p) assumes that at least some 
of % are nonzero. This phase transition is accompanied by the inversion symmetry breaking. 

The properties of thermodynamic quantities in the critical region, where the relative tempera¬ 
ture 

t = ——^ ->• 0 (4.6) 

is small, are characterized by the critical indices describing the behavior of the specific heat, 

Cv - T"“ , 


(4.4) 

(4.5) 


order parameter 

and isothermic compressibility 


(<P}~ t ? , 


K T OC X 


The dependence of an external field on the order parameter, at the critical temperature, is of the 
type 

hc<\(cp)\ s (T = T c ) . 

The pair correlation function, at large distance r = |r|, behaves as 


with the correlation length 


exp (-r/£) 

“ fd—2+ri ( r “»• °°) > 


And the vertex at T c exhibits the behavior 

r(ifc)« l+ck a (T = T c ) . 

Not all seven critical indices, a, /3, y, 8 , rj, v, ft), are independent. There are the so called 
scaling relations [53], due to Griffith, 


a + /3(l + <5) = 2 


and Widom, 


y+j8(l-5) =0, 


(4.7) 

(4.8) 
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from which the Rushbrook relation 

a+ 2/3 + 7 = 2 (4.9) 

follows. Also, the hyperscaling relations are known: 

a = 2 — vd, P = (d — 2 + 77 ) j , 

7=(2-??)v, 8 = d d + _\ + 1] r] • ( 4 - 10 ) 

Thus, only three critical indices, say rj, v, ft), can be treated as independent, and all others can be 
expressed through them. 

The critical indices rj, v, co can be represented [53] as expansions in powers of the variable 
e = 4 — d, formally valid for asymptotically small £ — > 0. We extrapolate such £-expansions by 
means of the factor approximants ( |3.35| ) and set e = 1 corresponding to d = 3-dimensional space. 
The results for all critical indices and different N are presented in the Table. The found values are in 
perfect agreement with experimental results, when these are available, and with numerical Monte 
Carlo calculations, as well as with complicated Pade-Borel summations, as discussed in [54]. It is 
interesting that for the limits N = —2 and N —> °o our method provides exact known results. 


Table: Critical indices for A-component <p 4 field theory. 


A 

a 

j 6 

7 

8 

V 

V 

CO 

-2 

0.5 

0.25 

1 

5 

0 

0.5 

0.80118 

-1 

0.36844 

0.27721 

1.07713 

4.88558 

0.019441 

0.54385 

0.79246 

0 

0.24005 

0.30204 

1.15587 

4.82691 

0.029706 

0.58665 

0.78832 

1 

0.11465 

0.32509 

1.23517 

4.79947 

0.034578 

0.62854 

0.78799 

2 

-0.00625 

0.34653 

1.31320 

4.78962 

0.036337 

0.66875 

0.78924 

3 

-0.12063 

0.36629 

1.38805 

4.78953 

0.036353 

0.70688 

0.79103 

4 

-0.22663 

0.38425 

1.45813 

4.79470 

0.035430 

0.74221 

0.79296 

5 

-0.32290 

0.40033 

1.52230 

4.80254 

0.034030 

0.77430 

0.79492 

6 

-0.40877 

0.41448 

1.57982 

4.81160 

0.032418 

0.80292 

0.79694 

7 

-0.48420 

0.42676 

1.63068 

4.82107 

0.030739 

0.82807 

0.79918 

8 

-0.54969 

0.43730 

1.67508 

4.83049 

0.029074 

0.84990 

0.80184 

9 

-0.60606 

0.44627 

1.71352 

4.83962 

0.027463 

0.86869 

0.80515 

10 

-0.65432 

0.45386 

1.74661 

4.84836 

0.025928 

0.88477 

0.80927 

50 

-0.98766 

0.50182 

1.98402 

4.95364 

0.007786 

0.99589 

0.93176 

100 

-0.89650 

0.48334 

1.92981 

4.99264 

0.001229 

0.96550 

0.97201 

1000 

-0.99843 

0.49933 

1.99662 

4.99859 

0.000235 

0.99843 

0.99807 

10000 

-0.99986 

0.49993 

1.99966 

4.99986 

0.000024 

0.99984 

0.99979 

OO 

-1 

0.5 

2 

5 

0 

1 

1 


13 
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5. Discussion 


We have considered the challenge of defining effective limits of divergent series by means of 
renormalization techniques. The necessity of this renormalization is dictated by the frequent occur¬ 
rence of such divergent series in complicated physical problems, e.g., arising when investigating 
phase transitions. The main ideas of optimized perturbation theory and self-similar approximation 
theory are presented. 


The name self-similarity comes from the self-similar relation (3.8) that is a necessary condition 


for the absolute minimum of the Cauchy cost functional (p.5|). The property of group self-similarity 
is equivalent to renormalization group in field theory [55]. To show this, let us introduce the 
variable 

(k = 0,1,2,...) . (5.1) 


T k = e k 


Instead of fl3.2| ), we can define 

y(Hk,<p) =Fk(x k ((p),Uk(x k ((p )) 

Since To = 1 at k = 0, the initial condition ( |3.3| ) becomes 

y{l,(p) = (p . 

In the place of the self-similar relation ([3T|), we now have 

= y{fk-,y{fpi ( p)) j 


(5.2) 


(5.3) 


(5.4) 


which defines a cascade. Embedding the cascade into a flow, according to ( p.l4[ ), with the discrete 
T k changing to the continuous t £ [0,°o), we obtain the scaling property 


y(MT,(p) = y(r,y(ju,<p)) . 


(5.5) 


This group property is typical of the renormalization group equations in field theory [56]. Instead 
of the Lie equation ( 3.16| ), we get 

dy(r ,(p) 


where the right-hand side 


dint 


m= 


= j8(y(T,<p)), 




(5.6) 


(5.7) 


T—\ 


is analogous to the Gell-Mann-Low function [57]. 

As a simple example of the application of self-similar approximation theory, let us briefly men¬ 
tion the calculation of the ground-state energy level for the one-dimensional anharmonic oscillator 
[58,59] with the Hamiltonian 

1 d 2 1 


H = ~2l? + 2 X ‘ + SX ' 


(5.8) 


where g £ [0,°°) is the anharmonicity, or coupling parameter. The introduction of a control function 
can be done through initial conditions, by starting perturbation theory with the Hamiltonian 

1 d 2 


„ IT 2 
H ° = ~2dx* + 2 X 


(5.9) 
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The control function Uk(g ) is defined by the quasi-fixed point condition (2.20). To first order 
of optimized perturbation theory, we have the energy E\ = E\(g). Using in the evolution integral 
( |3 .19 ), the cascade velocity v\ and t\ = 1 , we find for the self-similar approximation of the ground- 
state energy E* = E*(g) the equation 


4{E*) 2 — t _ e /_1_ J__±l 

4E?-1 P [4(E*) 2 -1 4E 2 -1 24 / 


(5.10) 


The comparison of E*(g ) with numerical calculations from the direct solution of the Schrodinger 
equation [60] shows that the found self-similar approximation E* (g) is applicable for all values of 

O.oo), yielding quite accurate results, whose maximal error does not exceed 0 . 1 %. 

In Sec. 4, we have illustrated the approach by calculating the critical indices for the inver¬ 
sion symmetry-breaking phase transition in an V-component rp 1 field theory. The results are in 
prefect agreement with experimental measurements as well as with complicated numerical tech¬ 
niques, such as Monte Carlo simulations or Pade-Borel summation. The advantage of our theory, 
as compared with other numerical methods, is the combination of much greater simplicity and very 
good accuracy. 
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